A  Yield  Strength  Model  and  Thoughts  on  an  Ignition 
Criterion  for  a  Reactive  PTFE-Aluminum  Composite 


by  M.  N.  Raftenberg,  M.  J.  Scheidler,  and  D.  A.  Casern 


ARL-RP-219 


August  2008 


A  reprint  from  the  Proceedings  of  the  55th  JANNAF  Conference ,  Newton,  MA,  12-16  May  2008. 


Approved  for  public  release;  distribution  is  unlimited. 


NOTICES 

Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position  unless 
so  designated  by  other  authorized  documents. 

Citation  of  manufacturer’s  or  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the 
use  thereof. 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the  originator. 


Army  Research  Laboratory 

Aberdeen  Proving  Ground,  MD  21005-5069 


ARL-RP-219 _ August  2008 


A  Yield  Strength  Model  and  Thoughts  on  an  Ignition 
Criterion  for  a  Reactive  PTFE-Aluminum  Composite 


M.  N.  Raftenberg,  M.  J.  Scheidler,  and  D.  A.  Casern 
Weapons  and  Materials  Research  Directorate,  ARL 


A  reprint  from  the  Proceedings  of  the  55th  JANNAF  Conference ,  Newton,  MA,  12-16  May  2008. 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering 
and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to 
comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


5c.  PROGRAM  ELEMENT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

A  reprint  from  the  Proceedings  of  the  55  th  JANNAF  Conference ,  Newton,  MA,  12-16  May  2008. 

14.  ABSTRACT 

We  studied  a  pressed  and  sintered  reactive  composite  of  74  wt%  polytetrafluoroethylene  (PTFE,  or  Teflon)  and 
26  wt%  aluminum  powder.  A  model,  which  we  call  “JCP”,  was  developed  to  relate  the  yield  strength  of  this  material  to  the 
equivalent  plastic  strain,  total  strain  rate,  and  temperature.  The  model  was  fit  to  Instron  compression  data  at  0.1/s  strain  rate 
and  split  Hopkinson  pressure  bar  (SHPB)  compression  data  at  approximately  2900/s  strain  rate.  The  SHPB  database  included 
initial  temperatures  of  297K  and  325K.  The  model  predicted  a  high  susceptibility  to  shear  localization  due  to  the  material’ s 
thermal  softening  and  strain  hardening  characteristics.  The  hypothesis  of  shear  localization  as  a  precursor  to  ignition  led  us  to 
consider  equivalent  plastic  strain  as  a  basis  for  an  ignition  initiation  criterion.  The  JCP  model  was  installed  into  the  EPIC  finite 
element  wavecode  and  used  to  simulate  Taylor  impact  tests  involving  impact  speeds  of  104  and  222  m/s.  The  computed 
boundary  shapes  versus  time  were  compared  with  digitized  images  from  a  framing  camera.  Agreement  was  good  at  1 04  m/s 
throughout  the  36-ps  time  range  of  observation.  Agreement  was  good  at  222  m/s  until  16  ps  after  impact.  Thereafter  the 
photographs  displayed  a  greater  degree  of  mushrooming  than  did  the  EPIC  simulation.  We  speculated  that  internal  ignition, 
fracture  of  various  modes  including  petalling,  and  inaccuracies  associated  with  interpolations  and  extrapolations  of  our  JCP  fit 
all  may  accounted  for  some  of  the  mushrooming  in  excess  of  that  predicted  using  JCP. _ 

15.  SUBJECT  TERMS 

reactive  material,  PTFE-aluminum,  Teflon,  ignition,  yield  strength 

17.  LIMITATION  18.  NUMBER 
OF  ABSTRACT  OF  PAGES 


UL  20 

Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


19a.  NAME  OF  RESPONSIBLE  PERSON 

M.  N.  Raftenberg 

19b.  TELEPHONE  NUMBER  ( Include  area  code) 

410-278-3684 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ARL-RP-219 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 


5d.  PROJECT  NUMBER 

622618H80 


1.  REPORT  DATE  (DD-MM-YYYY) 

2.  REPORT  TYPE 

August  2008 

Reprint 

4.  TITLE  AND  SUBTITLE 

A  Yield  Strength  Model  and  Thoughts  on  an  Ignition  Criterion  for  a  Reactive 
PTFE- Aluminum  Composite 

6.  AUTHOR(S) 

M.  N.  Raftenberg,  M.  J.  Scheidler,  and  D.  A.  Casern 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRD-ARL-WM-TD 
Aberdeen  Proving  Ground,  MD  21005-5069 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


A  YIELD  STRENGTH  MODEL  AND  THOUGHTS  ON  AN  IGNITION  CRITERION  FOR  A  REACTIVE 

PTFE-ALUMINUM  COMPOSITE 

M.  N.  Raftenberg,  M.  J.  Scheidler,  D.  A.  Casern 
U.S.  Army  Research  Laboratory,  Aberdeen  Proving  Ground,  MD 


ABSTRACT 

We  studied  a  pressed  and  sintered  reactive  composite  of  74  wt%  polytetrafluoroethylene  (PTFE, 
or  Teflon)  and  26  wt%  aluminum  powder.  A  model,  which  we  call  “JCP”,  was  developed  to  relate  the 
yield  strength  of  this  material  to  the  equivalent  plastic  strain,  total  strain  rate,  and  temperature.  The 
model  was  fit  to  Instron  compression  data  at  0.1/s  strain  rate  and  split  Hopkinson  pressure  bar  (SHPB) 
compression  data  at  approximately  2900/s  strain  rate.  The  SHPB  database  included  initial  temperatures 
of  297K  and  325K.  The  model  predicted  a  high  susceptibility  to  shear  localization  due  to  the  material’s 
thermal  softening  and  strain  hardening  characteristics.  The  hypothesis  of  shear  localization  as  a 
precursor  to  ignition  led  us  to  consider  equivalent  plastic  strain  as  a  basis  for  an  ignition  initiation  criterion. 
The  JCP  model  was  installed  into  the  EPIC  finite  element  wavecode  and  used  to  simulate  Taylor  impact 
tests  involving  impact  speeds  of  104  and  222  m/s.  The  computed  boundary  shapes  versus  time  were 
compared  with  digitized  images  from  a  framing  camera.  Agreement  was  good  at  104  m/s  throughout  the 
36-ps  time  range  of  observation.  Agreement  was  good  at  222  m/s  until  16  ps  after  impact.  Thereafter  the 
photographs  displayed  a  greater  degree  of  mushrooming  than  did  the  EPIC  simulation.  We  speculated 
that  internal  ignition,  fracture  of  various  modes  including  petalling,  and  inaccuracies  associated  with 
interpolations  and  extrapolations  of  our  JCP  fit  all  may  accounted  for  some  of  the  mushrooming  in  excess 
of  that  predicted  using  JCP. 


INTRODUCTION 

We  have  studied  a  pressed  and  sintered  stoichiometric  mixture  of  74  wt%  polytetrafluoroethylene 
(PTFE,  also  known  as  Teflon)  and  26  wt%  aluminum  powder.  Typical  average  pre-sintered  particle  sizes 
were  28  and  9  pm  for  the  PTFE  and  aluminum,  respectively.  This  material,  denoted  “PTFE-AI”,  has  been 
found  to  react  in  Taylor  impact  tests  at  sufficiently  high  impact  speeds.1,2  Figure  1  shows  impact  at  222 
m/s  of  a  rod  with  an  initial  length  of  50.8  mm  and  diameter  of  7.59  mm.  A  68-frame  Hadland  framing 
camera  was  used.  The  steel  anvil  appears  in  the  lower  portion  of  each  photograph,  and  an  aluminum 
sabot  appears  in  the  upper  portion.  This  sabot  had  a  mass  of  0.101  kg,  an  outer  diameter  of  39.6  mm, 
and  a  total  length  of  30.5  mm,  10.2  mm  of  which  was  hollowed  out  to  accommodate  the  rearmost  portion 
of  the  PTFE-AI  rod.  In  Figure  1  we  see  a  flash  near  the  mushroomed  leading  edge  of  the  specimen  at  32 
ps  post  impact. 

Figure  2  shows  a  conical  PTFE-AI  specimen  from  another  Taylor  impact  test.  This  test  involved 
an  impact  speed  of  104  m/s.  The  rod’s  initial  length  and  diameter  were  again  50.8  and  7.59  mm, 
respectively.  The  flat,  circular  base  of  this  cone  suggests  that  it  coincided  with  the  impact  face  of  the 
original  cylinder.  The  conical  surface  shown  in  Figure  2  displays  evidence  of  ductile  flow,  and  its 
blackened  regions  suggest  ignition.  On  the  basis  of  this  specimen,  we  developed  a  shear  localization 
hypothesis  and  adopted  a  metals-like  approach  to  modeling  mechanical  properties  of  PTFE-AI,  i.e.,  an 
equation  of  state  for  the  dilatational  response  and  von  Mises  plasticity  for  the  deviatoric  response. 

Reference  3  presented  data  from  uniaxial  compressive-stress  tests  on  PTFE-AI  with  an  Instron 
machine  at  a  constant  true-strain  rate  of  0.1/s  and  with  a  split  Hopkinson  pressure  bar  (SHPB)  at  an 
approximately  constant  true-strain  rate  of  about  2900/s.  The  Instron  specimen’s  initial  temperature  was 
296K.  Two  SHPB  specimens  were  initially  at  297K  and  two  were  at  325K.  These  Instron  and  SHPB  data 
are  collected  in  Figure  3. 

In  the  present  paper  we  propose  a  yield  strength  model  which  we  call  “JCP”,  an  explicit  function 
of  equivalent  plastic  strain,  equivalent-strain  rate,  and  temperature,  that  we  fit  to  the  data  in  Figure  3.  We 
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Figure  1.  A  Taylor  impact  test  shown  28  and  32  \is  after  initial  impact  at  222  m/s.  The  PTFE-AI 
cylinder’s  initial  length  and  diameter  were  50.8  and  7.59  mm,  respectively.  Note  the  steel  anvil  at 
the  bottom  and  the  aluminum  sabot  at  the  top  of  this  photograph.  (Reproduced  from  reference  3.) 
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Figure  2.  A  conical  PTFE-AI  specimen  recovered  from  a  Taylor  impact  test  at  104  m/s.  The  scale 
shown  is  in  millimeters.  The  cylinder’s  initial  length  and  diameter  were  50.8  and  7.59  mm, 
respectively.  (Reproduced  from  reference  3.) 
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Figure  3.  Instron  and  SHPB  Compression  Data  for  PTFE-AI.  (Reproduced  from  reference  3.) 


installed  this  model  into  the  EPIC  Lagrangian,  finite-element  wavecode4.  We  then  used  EPIC  to  simulate 
the  Taylor  impact  tests  at  104  and  222  m/s.  We  compared  computed  specimen  shapes  with  framing- 
camera  images  at  specific  times  post  impact.  At  the  higher  impact  speed,  we  obtained  general 
agreement  with  experiment  throughout  the  first  16  ps.  Thereafter,  JCP  under-predicted  the  extent  of 
specimen  mushrooming  near  the  anvil  interface.  We  attributed  this  discrepancy  to  one  or  more  of  three 
broad  causes:  interior  ignition,  material  damage,  and  deficiencies  in  our  representation  of  the  material’s 
undamaged,  inert  properties. 

In  the  remainder  of  the  paper,  we  provide  experimental  evidence  for  the  damage  mechanism 
known  as  “petalling”.  We  note  the  tendency  of  the  JCP  model,  with  its  explicit  thermal  softening,  to 
produce  the  phenomenon  of  adiabatic  shear  localization.  We  consider  the  latter’s  implications  for  PTFE- 
Al  ignition.  We  propose  a  model  that  can  introduce  some  effects  of  ignition  into  an  EPIC  simulation. 


RESULTS  AND  DISCUSSION 

THE  JCP  YIELD  STRENGTH  MODEL 


We  propose  the  following  functional  form  to  represent  the  constitutive  data  in  Figure  3. 


where 
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Here,  the  von  Mises  effective  stress,  <r,  is  an  explicit  function  of  equivalent  plastic  strain,  ep  ,  the  total 
equivalent-strain  rate,  s  ,  and  the  absolute  temperature,  9.  The  model  introduces  seven  material 
constants:  A0,  Au  B0,  B-i,  N0,  N-\  ,  and  9m. ,  plus  the  normalizing  strain  rate,  e0  . 


In  equation  1  thermal  softening  is  multiplicatively  decoupled  from  strain  hardening  and  strain-rate 
hardening.  The  material  constant  9m  is  the  parameter  governing  thermal  softening  and  is  the  temperature 
corresponding  to  zero  strength.  Note  that  equation  1  cannot  be  used  when  <9 exceeds  9m. 


In  equation  1,  A  and  C  are  two  independent  functions  of  total  equivalent-strain  rate  s.  The 
inverse  hyperbolic  sine  function  that  governs  the  strain-rate  dependence  is  related  to  the  logarithm 
function  by  the  identity 
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JCP  AND  A  CONSTANT-STRAIN-RATE.  CONSTANT-DENSITY.  ADIABATIC  PROCESS 

Suppose  a  JCP  material  undergoes  a  thermodynamic  process  that  is  adiabatic  and  involves 
constant  strain  rate  s, .  The  first  law  of  thermodynamics  can  be  used  to  derive  a  one-to-one  mapping 

between  temperature  and  equivalent  plastic  strain.  During  plastic  strain  increment  dsp ,  the  increment  in 
plastic  work  per  unit  mass  is 


dwP  -  <JJC?(£P  ,&,£*)  dsp 


Assume  that  any  work  associated  with  elastic  deformation,  whether  dilatational  or  elastic  deviatoric,  is 
relatively  negligible.  Assume  that  density  remains  approximately  constant  during  the  process,  so  that  the 
change  in  internal  energy  per  unit  mass  is 

de  -  p^cvd6 


where  cv  is  the  specific  heat  at  constant  volume.  Assume  that  the  Taylor-Quinney  factor  is  nearly  1 ,  so 
that  nearly  all  plastic  work  is  converted  into  heat.  The  first  law  of  thermodynamics  then  requires  that 

pQcvdO  -  crJCP  (sp dsp  (4) 


Combining  equations  (1)  and  (4), 
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Here  6}  is  the  initial  temperature.  The  result  is 


The  result  for  stress  in  the  material  governed  by  equation  1  during  an  adiabatic,  constant-density, 
constant-strain-rate  process  is 


where 
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FITTING  JCP  TO  THE  SHPB  AND  INSTRON  DATA 


The  total  strains  in  Figure  3  were  converted  to  plastic  strains,  spxx ,  by  means  of 
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P 

XX 


= 


(8) 


E  is  Young’s  modulus.  We  estimated  directly  from  the  Instron  curve  in  Figure  3  the  value  of  600  MPa  for 
E  at  low  strain  rate.  The  criterion  adopted  was  a  3%  strain  offset.  No  quantitive  estimate  for  E  could  be 
made  from  the  SHPB  loading  curves  because  equilibrium  was  not  attained  early  in  the  test.  We  obtained 
the  high-strain-rate  estimate  of  1 .6  GPa  for  E  based  on  unloading  curves  at  the  end  of  other  SHPB  tests. 
We  applied  equation  8  using  600  MPa  for  the  Instron  data  and  1 .6  GPa  for  the  SHPB  data,  and  obtained 
the  stress-plastic  strain  data  curves  in  Figure  4. 

In  the  case  of  uniaxial  stress  axx  with  all  other  components  zero,  the  von  Mises  stress  a  reduces 
to  <txx.5  If  incompressibility  is  assumed,  an  increment  d£p  of  equivalent  plastic  strain  reduces  to  d£px , 

the  plastic  strain  increment  along  the  direction  of  uniaxial  stress.5  A  similar  argument  can  be  used  to 
reduce  the  equivalent  total  strain  rate  to  the  component  along  the  direction  of  uniaxial  stress. 


The  four  SHPB  tests  in  Figure  3  all  involved  an  approximately  constant  strain  rate  of  2900/s.  We 
treated  these  tests  as  constant-volume,  adiabatic  and  applied  equations  6  and  7  with  s,  =  2900/s . 

Initial  temperature  9,  was  either  297K  or  325K.  We  obtained  from  reference  6  estimates  of  2270  kg/m3  for 
initial  density  p0  and  1161  J/kg-K  for  specific  heat  cv.  We  then  applied  the  FindFit  regression  function  in 

Mathematics7  software  to  fit  the  constants  A(2900/s),  B(2900/s),  N(2900/s) ,  and  9m.  The  results  are 
contained  in  Table  1 ,  and  the  JCP  fit  has  been  added  to  Figure  4. 

The  Instron  data  were  obtained  at  a  constant  strain  rate  of  0.1/s.  We  again  assumed  constant 
volume  and  adiabatic  conditions.  (The  latter  assumption  is  less  applicable  here  than  it  was  for  the  SHPB 
data.)  The  value  of  425K  obtained  for  9m  was  carried  over  from  the  SHPB  data  analysis.  Parameters 

,4(0.1/,s'),  5(0. l/s) ,  and  jV(0.  1/^)  were  evaluated  by  regression  using  FindFit  of  Mathermatica7 .  The 
results  are  given  in  Table  2.  The  Instron  fit  has  been  added  to  Figure  4. 

We  chose  0.03/s  for  the  normalizing  strain  rate  s0  because  it  is  substantially  smaller  than  the 

Instron  strain  rate  of  0.1/s.  Equation  2a,  evaluated  at  strain  rates  of  0.1  and  2900/s,  led  to  two  linear 
algebraic  equations  for  A0  and  A Similarly,  equation  2c  produced  two  equations  for  N0  and  /Vi.  B0  and 
Bt  were  then  evaluated  directly  from  equation  2b.  The  resulting  parameter  values  are  contained  in  Table 
3.  The  JCP  yield  strength  could  then  be  interpolated  and  extrapolated  in  terms  of  strain  rate  from  its  fits 
at  2900/s  and  0.1/s. 
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Figure  4.  JCP  Fits  to  SHPB  and  Instron  Compression  Data 


Table  1.  PTFE-AI  Material  Parameters  Used  in  or  Obtained  From  Fitting  JCP  to  SHPB  data 


P°  3 

(kg/m3) 

cv 

(J/kg-K) 

E 

(GPa) 

A( 2900 /s) 

(MPa) 

5(2900/s) 

(GPa) 

/V(2900/s) 

9m 

(K) 

2270 

1161 

1.6 

43.71 

1.998 

1.901 

425.3 

Table  2.  PTFE-AI  Material  Constants  Used  in  or  Obtained  From  Fitting  JCP  to  Instron  Data 


P°  3 

(kg/m3) 

cv 

(J/kg-K) 

E 

(MPa) 

i(0.1/s) 

(MPa) 

i(o.i/s) 

(MPa) 

jV(0.1/s) 

dm 

(K) 

2270 

1161 

600 

20.83 

57.34 

0.7718 

425.3 

Table  3.  JCP  Material  Constants  for  PTFE-AI 


(MPa) 

Ai 

(MPa) 

Bo 

(MPa) 

B-i 

No 

A/, 

4 

(1/s) 

dm 

(K) 

16.61 

2.227 

5.063 

23.21 

0.5633 

0.1099 

0.03 

425.3 

APPLICATION  OF  JCP  TO  THE  SIMULATION  OF  TWO  TAYLOR  IMPACT  TESTS 


The  JCP  model  was  installed  in  the  EPIC  Lagrangian  finite-element  wavecode4  and  applied  in 
simulations  of  two  Taylor  impact  tests  originally  reported  in  reference  1 .  In  reference  3  we  previously 
simulated  these  tests  using  the  Johnson-Cook  yield  strength  model8,  and  we  here  re-visited  them 
applying  JCP.  Both  rods  had  initial  length  and  radius  of  50.8  and  7.59  mm,  respectively.  They  were 
modeled  using  axisymmetric,  4-node  quadrilaterals  with  single-point  integration  and  an  initial  element  size 
of  0.1  mm.  The  aluminum  sabot  was  represented  by  concentrated  masses  distributed  along  the 
boundary  of  the  rearmost  1 0  mm  of  the  rod.  The  steel  anvil  was  modeled  as  rigid,  and  its  interface  with 
the  rod  as  frictionless.  The  dilatational  behavior  of  PTFE-AI  was  modeled  with  the  Mie-Gruneisen 
equation  of  state  representation  adopted  in  reference  6.  The  shock  Hugoniot  was  represented  by  a  linear 
shock  speed-particle  speed  relationship.  In  reference  6  an  intercept  of  1450  m/s  and  slope  of  2.258  were 
obtained  directly  from  shock-Hugoniot  data.  The  constant  Gruneisen  parameter  was  estimated  to  be  0.9. 
The  elastic  shear  modulus  was  554  MPa,  which  we  determined  from  the  high-rate  Young’s  modulus  of 
1.6  GPa,  and  the  elastic  bulk  modulus  consistent  with  the  1450-m/s  intercept  of  the  shock  Hugoniot. 

Figure  5  shows  the  computed  rod  shapes  for  104-m/s  impact  superimposed  on  the  boundary 
locations  digitized  from  framing-camera  photographs.  Note  that  radial  displacements  have  been 
amplified  by  a  factor  of  5  in  this  figure.  The  EPIC  simulations  agree  reasonably  well  in  terms  of  the  rod’s 
shape  versus  time.  The  final  outer  radius  at  the  anvil  interface  is  in  particularly  close  agreement. 

Figure  6  shows  computed  contours  of  equivalent  plastic  strain  and  temperature  at  36  ps  after 
impact  at  1 04  m/s.  Note  the  conical  shape  of  these  contours.  If  we  adopt  the  hypothesis  that  an  ignition 

criterion  should  be  based  on  equivalent  plastic  strain  (e^  ),  then  the  conical  contours  in  Figure  6  support 
our  earlier  speculation  that  ignition  first  occurred  on  the  surface  of  the  recovered  specimen  in  Figure  2. 

Figure  7  compares  with  experiment  the  computed  specimen  shape  for  various  times  after  222- 
m/s  impact.  (Radial  displacements  are  not  multiplied  by  5  as  they  are  in  Figure  5.)  The  computed  shape 
shows  good  agreement  with  experiment  at  4,  8,  12,  and  16  ps  after  impact.  By  20  ps,  however,  the 
mushroomed  volume  is  substantially  larger  in  the  experiment  than  in  the  calculation.  This  discrepancy 
between  theory  and  experiment  becomes  successively  more  pronounced  at  24,  28,  and  32  ps. 

Thus,  at  an  impact  speed  of  222  m/s,  the  EPIC  simulation  under-predicts  the  extent  of 
mushrooming  at  20,  24,  28,  and  32  ps  but  not  at  4,  8,  12,  and  16  ps.  Several  explanations  can  be 
envisioned.  First,  the  JCP  fit  might  not  apply  throughout  the  ranges  of  strain,  strain  rate,  and  temperature 
that  the  Taylor  rod  underwent  in  the  course  of  the  test.  Specifically,  the  fit  was  interpolated  in  terms  of 
strain  rate  and  temperature  and  extrapolated  in  terms  of  strain,  strain  rate  and  temperature.  Second, 
there  may  have  been  fracture  and  associated  void  formation  within  the  rod.  Third,  there  may  have  been 
interior  ignition  and  associated  gas  formation. 


DIRECT  OBSERVATIONS  OF  PETALLING 


Figure  8  contains  photographs  from  a  modified  Taylor  impact  test,  using  a  small-diameter  steel 
anvil.  The  PTFE-AI  rod’s  initial  length  and  diameter  were  25.4  mm  and  6.4  mm,  respectively.  The  anvil’s 
diameter  was  9.5  mm.  The  impact  speed  was  202  m/s.  The  photograph  at  50  ps  clearly  displays  the 
failure  mode  known  as  “petalling”,  which  refers  in  this  case  to  radial  cracks  in  the  mushroomed  region. 
Retailing  involves  tensile  failure.  In  the  mushroomed  region,  the  rod  undergoes  progressive 
circumferential  tension,  eventually  leading  to  rupture. 

CRITICAL  STRAIN  FOR  SHEAR  LOCALIZATION:  IMPLICATIONS  FOR  IGNITION 

In  Figure  10,  we  have  extrapolated  to  larger  strain  values  the  JCP  fit  in  Table  1 .  The  JCP  fit, 
evaluated  for  either  initial  temperature,  attains  its  maximum  value  at  an  equivalent  plastic  strain,  sp  ,  of 

0.66.  As  strain  is  increased  above  sp  ,  the  decrease  in  due  to  thermal  softening  has  a  larger 

magnitude  than  the  increase  due  to  strain  hardening.  The  resulting  negative  slope  in  the  adiabatic  stress- 
strain  curve  at  strains  greater  than  sp  can  lead  to  the  phenomenon  of  shear  localization  such  as  occurs 

in  certain  metals.  The  material  is  locally  unable  to  equilibrate  its  load.  Additional  plastic  flow  generates 
an  additional  temperature  rise,  which  further  softens  the  material.  A  self-sustaining  process  develops 
which  leads  to  the  local  state  of  zero  deviatoric  stress  and  the  plateau  value  of  6m  for  temperature. 
Temperature  as  a  function  of  strain  is  given  by  equation  5  and  has  been  added  to  Figure  10. 

The  shear  localization  phenomenon  may  have  an  added  significance  in  the  case  of  PTFE-AI,  with 
its  potential  for  ignition.  We  speculate  that  the  concentrated  shearing  motion  may  be  able  to  strip  the 
oxide  layer  from  aluminum  particles  in  the  shear  band,  thereby  exposing  those  particles  to  the 
surrounding  PTFE.  Perhaps  such  exposure  is  a  necessary  condition  for  ignition. 


The  critical  strain  sp  is  defined  by 


d<?jcp>  (,£P  ,£•) 


dsp 


=  0 


Substitution  of  equation  6  into  9  yields  the  algebraic  equation 


(£-p)2"  i  294K)  (,^y_i  + 

B  B 


=  0 


(9) 


(10) 


The  coefficients  of  this  equation  involve  A  ,  B  ,  and  N ,  which  are  functions  of  strain  rate.  It  follows  that 
the  solution  for  i:’’  is  also  a  function  of  strain  rate.  Hence,  if  ignition  is  initiated  by  shear  localization,  we 
expect  that  the  value  of  strain  corresponding  to  ignition  initiation,  call  it  e?  ,  is  a  function  of  strain  rate. 
Moreover,  the  physical  processes  that  intervene  between  shear  localization  and  ignition  can  introduce 
additional  time  scales  and  thereby  make  more  strongly  dependent  on  strain  rate  than  is  ep  .  These 

processes  include  the  mechanical  stripping  of  the  aluminum  particles’  oxide  layer  and  startup  of  the 
chemical  reaction(s). 

A  MODEL  FOR  PARTIAL  EFFECTS  OF  IGNITION 


A  model  that  introduces  into  a  PTFE-AI  impact  simulation  some  effects  of  ignition  was  added  to 
EPIC  to  use  in  conjunction  with  JCP  yield  strength.  This  ignition  effects  model  introduces  three  material 

constants:  ignition  strain,  epn ,  potential  chemical  energy  per  unit  mass,  e0 ,  and  the  time  duration  of  the 
chemical  reaction,  Arign .  The  model  consists  of  two  parts:  a  criterion  for  detecting  ignition  and  a  post- 
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Figure  5.  Experimental  and  computed  specimen  boundaries  in  the  Talyor  impact  test  at  an  impact 
speed  of  104  m/s.  (Experimental  data  reproduced  from  reference  3.) 


Figure  6.  Contours  of  equivalent  plastic  strain  and  temperature  at  36  ps  after  impact  at  104  m/s. 
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Figure  7a.  Experimental  and  computed  specimen  boundaries  in  the  Taylor  impact  test  at  an 
impact  speed  of  222  m/s  (4,  8, 12,  and  16  ps  after  impact).  (Experimental  data  reproduced  from 
reference  3.) 
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Figure  7b.  Experimental  and  computed  specimen  boundaries  in  the  Taylor  impact  test  at  an 
impact  speed  of  222  m/s  (20,  24,  28,  and  32  ps  after  impact).  (Experimental  data  reproduced  from 
reference  3.) 


Figure  8.  A  modified  Taylor  impact  test  with  an  PTFE-AI  specimen  impacting  at  202  m/s  a  small- 
diameter  steel  anvil.  Retailing  fracture  has  begun  by  20  ps  and  is  pronounced  at  52  ps. 
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Figure  10.  The  JCP  Fit  ats  =  29001  s ,  Extrapolated  to  Larger  Values  of  Strain.  Note  the  Critical 
Strain  of  0.66  Corresponding  to  Maximum  Stress. 

detection  algorithm.  The  detection  criterion,  applied  individually  to  each  finite  element,  is 


>£p 
—  ign 


(11) 


The  element-dependent  parameter  rign  is  the  time  associated  with  the  time  step  at  which  equation  1 1  was 

first  satisfied  within  that  element.  The  element’s  increment  in  specific  internal  energy  (per  unit  mass)  is 
computed  for  a  given  time  step  of  duration  At  at  time  t  using 


Aeign  = 


At 

A  t-„„ 


t<t:r 


;  tisa<t<tiw+Atk 


’  ^ign  A^ign  ^  t 


(12) 


Aeign  is  identified  with  the  increment  in  specific  internal  energy  generated  by  the  chemical  reaction. 

The  post-detection  algorithm  applies  Aeign  as  a  source  term  in  the  first  law  of  thermodynamics, 


PqCvA@  <TjCP(£  ,  0 *  £)  Ac  +  Aeign 


and  in  the  Mie-Griineisen  equation  of  state, 


P(P,e) 


(1-S?7)2l  2] 


+  PPo  (e  +  Aeign) 


(13) 


(14) 


In  equation  1 1  A# and  A sp  are  the  increments  in  temperature  and  equivalent  plastic  strain,  respectively, 
in  that  element  during  the  given  time  step.  In  equation  12,  ij  =  1  -  /?0/ p ,  where  p  is  the  current  density 
within  the  element. 


Equation  1 1  contains  the  assumptions  that  the  process  is  effectively  adiabatic  and  that  the  energy 
increment  associated  with  the  material’s  density  change  is  negligible  relative  to  the  increments  in  plastic 
work  and  chemical  energy.  According  to  equation  1 1 ,  the  introduced  chemical  energy  raises  the 
element’s  temperature.  This  reduces  the  element’s  yield  strength  by  equation  1 .  In  a  Taylor  impact 
simulation,  a  result  would  tend  to  be  an  increased  mushrooming.  According  to  equation  12,  the 
introduced  chemical  energy  is  accompanied  by  a  rise  in  pressure.  This  would  produce  an  increased 
specific  volume  and  also  tend  to  lead  to  an  increased  mushrooming. 

The  model  ignores  a  known  effect  of  ignition,  namely  the  material’s  phase  change  into  gaseous 
products.  At  a  given  time  following  the  local  onset  of  ignition,  the  material  within  a  lagrangian  finite 
element  should  be  divided  between  a  certain  volume  fraction  still  in  the  solid  phase  and  the  remaining 
volume  fraction  converted  to  a  gas.  Equations  13  and  14  should  properly  only  apply  to  the  material  still  in 
the  solid  phase. 


SUMMARY  AND  CONCLUSIONS 

The  JCP  model  was  developed  to  relate  the  von  Mises  yield  strength  of  a  reactive  PTFE- 
aluminum  composite  to  equivalent  plastic  strain,  total  strain  rate,  and  temperature.  The  model  was  fit  to 
Instron  and  SHPB  compression  data  at  strain  rates  of  0.1/s  and  approximately  2900/s,  respectively.  The 
model  was  installed  into  the  EPIC  finite  element  wavecode  and  used  to  simulate  Taylor  impact  tests  at 
impact  speeds  of  104  and  222  m/s.  The  rods’  initial  length  and  diameter  were  50.8  and  7.59  mm, 
respectively.  The  computed  boundary  shapes  versus  time  were  compared  with  digitized  images  from  a 
framing  camera.  Agreement  was  good  at  104  m/s  throughout  the  observed  36-ps  duration.  At  222  m/s, 
agreement  was  good  until  16  ps  after  impact.  Thereafter  the  photographs  displayed  a  greater  degree  of 
mushrooming  than  did  the  EPIC  simulation.  Our  speculated  reasons  for  this  under-predicted  late-time 
mushrooming  included  inaccuracies  introduced  by  extrapolations  and/or  interpolations  of  our  JCP  fit, 
material  fracture,  and  ignition.  The  fracture  hypothesis  was  then  supported  in  a  modified  Taylor  impact 
test  in  which  petalling  fracture  was  observed  by  means  of  a  higher  resolution  camera. 

A  model  was  developed  to  introduce  partial  effects  of  ignition  into  an  EPIC  simulation.  The 
ignition  criterion  is  based  on  a  critical  value  of  equivalent  plastic  strain.  After  detection,  effects  of  the 
chemical  reaction  are  introduced  by  means  of  a  source  term  in  the  first  law  of  thermodynamics  and  in  the 
Mie-Griineisen  equation  of  state.  The  model  is  partial  in  that  transformation  into  gaseous  products  has 
not  yet  been  taken  into  account.  This  model  should  be  refined  and  then  applied  in  Taylor  impact 
simulations.  A  model  for  petalling  fracture  is  also  required  for  a  more  complete  Taylor  impact  simulation. 
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RESEARCH  TRIANGLE  PARK 
NC  27709-2211 


NO.  OF 

COPIES  ORGANIZATION 


ABERDEEN  PROVING  GROUND 


42  DIR  USARL 

AMSRD  ARL  WM 
R  DOWDING 
J  MCCAULEY 
S  MCKNIGHT 
J  SMITH 

AMSRD  ARL  WM  B 
J  NEWILL 

AMSRD  ARL  WM  BD 
B  FORCH 
B  HOMAN 
K  MCNESBY 
AMSRD  ARL  WM  MB 
D  HOPKINS 
L  KECSKES 
AMSRD  ARL  WM  MD 
B  CHEESEMAN 
ECHIN 
G  GAZONAS 
C  RANDOW 
C  YEN 

AMSRD  ARL  WM  T 
P  BAKER 

AMSRD  ARL  WM  TB 
R  BANTON 
R  BITTING 
W  BUKOWSKI 
N  ELDREDGE 
R  GUPTA 
S  KUKUCK 
AMSRD  ARL  WM  TC 
T  FARRAND 
L  MAGNESS 
D  SCHEFFLER 
B  SCHUSTER 
S  SEGLETES 
AMSRD  ARL  WM  TD 
S  BILYK 
T  BJERKE 
D  CASEM 
J  CLAYTON 
Y  HUANG 
R  KRAFT 
BLOVE 

M  RAFTENBERG  (5  CPS) 
M  SCHEIDLER  (2  CPS) 

T  WEERASOORIYA 


